Constraints on the ubiquity of coronal X-ray cycles 
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ABSTRACT 

Stellar activity cycles are known to be a widespread phenomenon amongst moderately active 
solar- and late-type stars from long-term periodic variations in chromospheric Ca II H and K 
emission lines, yet to date only a handful of coronal X-ray cycles are known. We have surveyed 
serendipitously observed stellar sources in fields observed multiple times in the last decade by 
XMM-Newton and present our analysis of 9 stars from 6 fields. Since our sample is flux-limited, 
it is strongly biased towards higher levels of X-ray activity. We fit a single temperature APEC 
spectrum to each source and search for significant periodicities using a Lomb-Scargle Periodogram 
(LSP). We use a Monte Carlo (MC) algorithm to yield robust analysis of the statistical significance 
of cycle detections and non-detections. None of the 9 stellar lightcurves show any convincing 
indications of periodicity. From MC simulations we simulate the detection capabilities of our 
methodology and, assuming a uniform distribution of cycle periods and strengths over the domain 
searched, we conclude with 95% confidence that less than 72% of the stars represented by our 
sample of active stars have 5-f 3 year coronal X-ray cycles. 

Subject headings: Stars: activity, stars: coronae, X-rays: stars 



1. Introduction 

The 1 1 year s un-spot cycle wa s discovered over 
160 years ago bv ISchwabe ( 1843 ). Later, the sub- 
sequent discovery by Halej (|l908l ) that the global 
solar magnetic field reverses polarity after each 
sun-spot cycle led to our current understanding 
of the 11 year solar activity cycle as arising out of 
an underlying 22 year magnetic cycle. 

In 1966, in order to ascertain whether or not 
other stars exhibited long-term periodic behav- 
ior, O. C. Wilson started a long-term monitor- 
ing program of chromospheric activity in solar- 



type main sequence stars. He defined the S-index, 
which characterizes the strength of the Ca II H 
and K lines re lative to the photospheric continuum 
(|Wilsonlll9 68). Through monitoring a wide range 
of stars in the chromospheric Ca II H and K emis- 
sion lines via the S-index, we now understand that 
the activity cycle of our Sun is a mor e general phe- 
nomenon amongst s olar type stars ( Wilson! 1978t 



Baliunas et al.lll995l ) . What factors determine the 



1 Based on observations obtained with XMM-Newton, 
an ESA science mission with instruments and contributions 
directly funded by ESA Member States and NASA. 



durations of these cycles, however, remains elusive 
to theoretical explanation. 

The Ca n H and K lines are well established 
indicators of magnetic activity. Additionally, vir- 
tually all solar- and late-type main sequence stars, 
including the Su n, have hot coronae that emit in 
the X-ray bands (jVaiana et al.lll98ll ). The energy 
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required to heat the corona is supplied by the stel- 
lar magnetic field, the dynamics of which are gov- 
erned by an a — Q dynamo involving differential 
rotat ion and turbu l ence in the convectiye env elope 
(e.g. iRosner et al.l Il985t lOssendriiver 20031 and 
references therein), which is clearly manifested in 
the stellar activity — rotation relationship: stars 
that rotate faster are typically more X-ray lu- 



mino us (e.g.. iPallavicini et al.lll981t IWright et al 



2011). 



Early efforts to statistically determine a trend 
of long-term variability in X-ray stars from 
Einst ei n and ROSAT data are summarized in 
ISternl ( 19981 ). Stars that are X-ray faint typi- 
call y show low levels of short-timescale var iabil- 
ity (|Feigelson et al.ll2004t IWright et aLll2010h. but 



may v ary on longer timescales. Micela fc Marinol 
((2003J) compare the solar variability to stellar data 
and find that this is consistent with cycles ex- 
plaining much of the apparent luminosity spread 
in solar-liked stars. Stars that are X-ray bright 
(Lx > 10 28 erg s _1 ) are not believed to exhibit 
long term X-ray cycles and instead have X-ray lu- 
minosities that vary erratic ally on short timescales 
(see review by IStern 19981 ). In this paper, we add 
experimental credence to this assertion by demon- 
strating an absence of long-term cycles in X-ray 
active solar-like stars. 

Currently, long term X-ray cycles have been 
found in only a handful of stellar sources: 61 
Cvgni A dHempelmann et al.l 120061). HP 81809 
(iFavata et al.ll2008l). a Centauri B ([Robrade et al 



2005: lAvres et al.ll2008t |A^rej|2009|) and our Sun 



The amplitudes of the stellar X-ray cycles in these 
objects are typically larger than their respective 
S-index cycles by up t o an order of magn itude, as 
is the case in the Sun ([Judge et al. I l2003h . 

The large levels of periodic X-ray variability ex- 
pected from solar-like main-sequence stars mean 
that their cycles should be detectable with mod- 
ern X-ray observatories such as XMM-Newton that 
have been in service for over a decade. In this ar- 
ticle we analyze serendipitous stellar sources de- 
tected in XMM-Newton fields that have been fre- 
quently observed over the lifespan of the instru- 
ment (~12 years) using publicly available data 
from the XMM-Newton archive. In Section [5] we 



discuss the observations used for this study and 
the methods employed. In Section |3] we present 
the observational results and Monte Carlo (MC) 
simulations to quantify our detection limits. In 
Section 2] we discuss the implications of these re- 
sults on our understanding of stellar X-ray emis- 
sion and long-term, solar-like X-ray cycles in other 
stars. We end with a presentation of our conclu- 
sions in Section [5] 

2. Observations and Data Analysis 

In XMM-Newton's nearly 12 year history, nu- 
merous fields have been observed semi-regularly 
for purposes ranging from calibration (e.g. Capella) 
to measuring long term periodic variability in X- 
ray luminosities (e.g. HD 81809). These fields 
offer a rare opportunity for the detection of long 
term stellar X-ray cycles, which are believed to 
correlate strongly with chromospheric activity cy- 
cles. While centr al sources have be en analyzed 
(e.g., HD 81809. iFavata et al.l [200ll . serendipi- 



2 "solar-like" refers to spectral type (F-M), and not levels of 
X-ray activity. 



tous sources have not, to our knowledge, been 
checked for periodic behavior. Thus, for fields 
which have a long and consistent history of XMM- 
Newton observations, we check for periodicities in 
the lightcurves of serendipitous stellar sources. 

2.1. Field selection 

Target fields that have been frequently observed 
(i.e. that have a baseline of at least six years and 
have been observed during at least 5 semesters be- 
tween 1999 and 2012) were selected for analysis. 
We will show below that all candidate sources are 
relatively faint. To provide enough counts for a 
stable fit of the spectrum, we only use observa- 
tions with a Good Time Interval (GTI) of >10 ks 
for one or more of the MOS1, MOS2, and PN cam- 
eras. Any data mode is acceptable. For instance, 
even though modes PPW2 and PPW4 have small cen- 
tral windows that do not read out parts of the field 
near the central source, these exposures still con- 
tain most serendipitous sources in the field. Six 
fields pass these selection criteria. Table [1] sum- 
marizes the observation timing and baselines for 
the 6 fields used in this paper. "Distinct obser- 
vations" are the number of half-year intervals be- 
tween 1999 and 2012 containing at least one usable 
observation. 
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2.2. Source selection and identification of 
stars 

X-ray sources are selected by cross- ma tching 



the 2 XMMi Serendipitous Source Catalo gue ([W atson ct 
2009) with the optical NOMAD catalog (jZacharias et al 
2004 ), using a matching radius of 5"consistent 
with the XMM-Newton PSF for off-axis sources. 
For bright central sources, the 2XMMi catalog 
contains many spurious detections in the vicinity 
that are related to read-out streaks in the orig- 
inal data. These read-out streaks are manually 
removed from our target list. 

We impose a minimum flux limit of 25 counts 
in a 10 ks observation, which equates to a maxi- 
mum flux uncertainty of 20% from Poisson statis- 
tics. For an observation of 10 ks and 25 counts this 
equates to a flux limit of 1.4 x 10 -14 erg cm -2 s _1 
for a typical stellar observation and stellar prop- 
erties (solar abundance, plasma temperature of 
1 keV and observations using the medium filter 
for the MOS). 

A number of steps are taken to identify stellar 
X-ray sources. First, because stars have relatively 
soft X-ray spectra compared to most extragalactic 
sour ces, we require a ha rdness ratio of HR1 > 
(see Watson et al. 2009t ). We then search the lit- 
erature for available photometry and spectroscopy 
to allow us to classify sources as stellar or non- 
stellar. Due to the large number of non-stellar 
sources in deep X-ray observations, we require a 
high-confidence verification of the stellar nature 
of our candidates, discarding all sources that lack 
sufficient data to confidently classify them. We 
use three methods for identifying stellar sources 
based on the availability of certain observations 
and measurements, which, in order of priority, are: 



Table 1: Field Characteristics 



al 



• Spectra showing that the source is stellar 
in nature. Available for five targets (see Ta- 
ble E). 

• Proper motions showing that the source is 
sufficiently nearby as to have a measurable 
motion on the plane of the sky with a signif- 



Central 


Abbr. 


Baseline 


Distinct 


source 


name 


(yr) 


obs. 


Capella 


CAP 


7.48 


6 


Zeta Puppis 


ZP 


9.98 


13 


PKS 2155-304 


PKS 


10.9 


16 


3C273 


3C 


10.49 


11 


HD 81809 


HD 


6.0 


5 


Mkn 421 


MK 


10.98 
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icance of 3a from the catalog of iRoser et al 
(|2008h . For stars with proper motion, but 



no spectrum, we require that the photom- 
etry -if available - is also consistent with 



stellar- like colors ([Covey et al.ll2007n . Three 
targets meet these criteria (#4, #5, #7). 

• Photometry indi cating stellar-like colors 
(jCovev et al.l 120071 ) in all bands. One tar- 
get meets this criterion (#3). 

In addition to these requirements, stellar tar- 
gets without available spectral information are re- 
quired not to have been identified as extragalac- 
tic by any other studies incorporating observations 
at other wavelengths, morphology, or more com- 
plex multi- wavelength diagnostics. These crite- 
ria reduce the initial list of 72 sources to a high- 
confidence list of 9 stellar sources (see Table [2]) . 

2.3. Spectral types and distances 

For these 9 stellar sources, spectral types are 
determined from the best available data, either via 
spectroscopy (for 5 sources) or from the photomet- 
ric color with the largest baseline ( for 4 sources) 
and u sing the empirical data of IPecaut et al 
(|2012[ ) . All sources are assumed to lie on the main- 



sequence, a reasonable assumption since our tar- 
geted fields are away from known star forming re- 
gions and X-ray emission from la te-type post-main 
seque nce stars is rare and weak (|Linskv fc Haischl 
19791) . The determined spectral types range from 
F-type to early M-type. Extinctions were calcu- 
lated based either on the available photometry 
(for stars with known spectral types), or from 
an initial estimate of the distance (from the red- 
dened magnitudes) and assuming an extinction of 
Ay — 1 mag / kpc, up to a maximum extinction 
of the total integ rated Galactic extinc tion in that 
line-of-sight from lSchlegel et al.l (119981). Distances 



were determined using 2MASS (jSkrutskie et al 
20061 ) if s -band photom etry and the empiric al ab- 
solute magnitudes from IPecaut et al.l (|2012t ) . For 
one source, Gliese 195 A, a parallax is available 
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(|Jenkinslll952f) and is therefore used to determine 
its distance. The determined stellar distances 
range from 14 to 663 pc. 

2.4. Data analysis 

Observational data are obtained from the 
XMM-Newton archive and analyze d using SAS 
(version 11.0.0, iGabriel et all 120061 ). We apply 
all standard selection criteria for rejecting cos- 
mic rays and filtering out the background rate to 
avoid contamination from solar proton flares. For 
all sources, spectra are extracted and the appro- 
priate response files are generated. Fitting of spec 



tral m odels is performed in Sherpa (jFreeman et al 
20 llh . 



Due to the rather large point spread function 
of XMM-Newton and the relatively dim nature of 
the selected sources, a 450-500 pixel (~8-9") ra- 
dius is used for the source region and a 2500 pixel 
(~46") radius is used for the background region, 
manually chosen to avoid other 2XMMi sources. 
One background region is chosen for each field to 
simplify the processing^ 

2.5. Fitting APEC models 

A single temperature APEC model is fitted 
to extracted spectra that contain a minimum 
of 125 counts within the 0.4-7 keV range af- 
ter subtracting background counts. We fix 
the interstellar absorption using the Ay val- 
ues from table [2] and the relation Nh = 1.8 X 



1 ° 21 mag j crfF A v dPredehl fc Schnhttl Il995h . We 
employ the Levenberg-Marq uardt (LM) fitting 



algorithm ( Marquardt 19631 called le vmar in 
sher pa) with a Gehrel's \ 2 statistic (jGehrels 
19861 ) to fit APEC models to the observed spec- 
tra. If the Gehrel's \ 2 value of a LM fit exceeds 
a threshold value heuristically set to 0.7, we seek 
an improvement of the fit. First, the LM algo- 
rithm is again tried with five arbitrary and dis- 
tinct initial parameter values. If the Gehrel's x 2 
of the fit is still above the heuristic threshold, 
the fitting method is changed to the more ac- 
curate but slower d ifferential evolution method 
( Storn &: Price 1997? ) (called moncar in sherpa) 
as a second effort for improving the fit. 



At the end of the fitting procedure, if the 
Gehrel's \ 2 lS below 0.25 or above 3, la flux uncer- 
tainties are not computed. A \ 2 above 3 strongly 
suggests that the model does not adequately fit 
the data. A x 2 below ^0.25 suggests that the 
APEC model fit is not well constrained. Flux values 
without estimated la uncertainties are not used in 
later analysis. 

Unless otherwise stated, X-rays are defined to 
be in the 0.4-7 keV band. Fx is the integrated 
energy flux over this range. Lx is the unabsorbed 
luminosity over this range. The Fx values are ob- 
tained from integrating the fitted APEC model. Fit 
errors for the integrated Fx are calculated from 
the uncertainty of the model normalization. 

3. Results 

3.1. The stellar sample 

Table summarizes the properties of our stel- 
lar sample. Uncertainties in Lx and Lx/Lboi are 
about 0.2 dex, including the statistical uncertainty 
from the fitted X-ray model and the uncertainty 
in distance from the photometry. The spectral 
types of our sample range from F to early M, the 
activity level (log Lx/Lb i) from —3 to —6. Ac- 
tive stars are known to saturate at an activity 
level of \og(Lx / Lboi) ~ —3 , while the activity 
level drops for older or more slowly rotating stars 
( Stauffer et al.lll994T ). Naturally, stars as inactive 
as our Sun are not part of this serendipitous sam- 
ple because they are too X-ray faint. 

3.2. The Lomb Scargle Periodogram 

After fitting a single temperature APEC model 
to the spectra of all 9 stellar sources and integrat- 
ing the fitted model over the 0.4-7 keV range to 
obtain flux estimates for each exposure, the re- 
sulting lightcurve is analyzed for the presence of 
periodic behavior. This is done using a Lomb- 
Scargle Periodogram (LSP) (iLombl 119761 : IScargld 
19821 ). Derived quantities from the LSP are (1) 



3 Using background regions from different chips yielded con- 
sistent results in a small test run. 



the frequency of the best fit sinusoidal function 
to the data (known as the peak frequency) and 
(2) the "False Alarm Probability" (FAP) of the 
strongest periodicity found. The FAP is defined 
to be the probability that periodic behavior of a 
certain strength could be found in Gaussian noise. 
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Table 2: Summary of stellar sources 
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(pc) 


(erg/s) 




1° 


CAP 
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46:08:27.6 


F0 d 


8.06 


7.23 
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2<> 


CAP 


05:17:24.0 
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ZP 
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GO 


13.37 
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HD 
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K3 
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HD 
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0.11 
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-3.9 
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MK 
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-4.4 
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3C 
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14.09 
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0.06 
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Other names: a BD+45 1076; fc Capella H, G 96-29, Gliese 195 A; C HD 95976 

Notes: Spectral types e xtracted fr o m the literature where a citation is provided, or from a vailable colors. V 
magnitudes taken from Hog et all d2000T) and the NOMAD catalog (jZacharias et al.ll2004f) . K s magnitudes 
extracted from the 2MASS catalog (jSkrutskie et al.l 120061 ) and distances determined photometrically from 
this (with the exce ption of #2 , w hich h as a p arallax) 

Re ferences: teeckmannl (Il975 ^ 

Lopez-Santiago et al.l (|2007l h ICovev fc et al 



] Mo ntes et al 
" (|2008h 



d200lh. - ICannon fc Pickering! 
Jenkins! (|l952h . 



mm. 



3.3. Cycle detection 

Bali unas et al. | (|l995f ) consider a FAP < 0.001 



(>99.9% confidence) as a detection of periodicity 
in their analysis of S-index data from the Mount 
Wilson Observatory. In our case, X-ray flux un- 
certainties are large and most stars have been ob- 
served sparsely in the last 12 years. Thus, we do 
not expect the lightcurves of the stellar sources in 
this paper to produce such high-confidence statis- 
tics. 

Instead of employing a simple FAP cutoff, we 
use a MC algorithm to assess the statistical like- 
lihood of LSP results by taking into consideration 
the flux uncertainties, which vary from observa- 
tions to observation. We consider a lightcurve to 
be periodic (or at least warrant further investiga- 
tion) if 68% of all MC-sampled lightcurves have a 
FAP below 0.30. 

3.4. Monte Carlo simulations 

While periodogram analysis offers an invaluable 
statistical method for determining the presence of 
periodicities, the LSP calculation itself does not 
account for uncertainties in the data. Therefore, 
we use a MC algorithm to determine the statistical 
significance of a detected periodicity. 

In a lightcurve containing T data points with 
uncertainties, the MC algorithm samples T flux 



values in the following way for each step (N = 
1000 total steps) of the simulation: 

FMc(ti) — Fi + Xi , 

where Fi is the ith observed flux (at time U) and 
Xi is a random variable sampled from a Gaus- 
sian distribution of mean and a standard devia- 
tion equal to the uncertainty of the empirical flux 
value. Figure Q] shows the lightcurve, a cumulative 
histogram of FAP values and a histogram of the 
strongest periods (inverse of the peak frequency in 
the LSP) detected in #7, a Kl star in the 3C 273 
field. 

3.5. The test case: HD 81809 



As a test case, we present an analysis of a stel- 
lar source with a known X- ray cycle. HD 81809 
has an 8.2 year S-index cycle (|Baliunas et al ll995h 
and an X-ray cycle that appears in phase with the 
S-index cycle (jFavata et al.ll20081 ). We use all pub- 
lically available XMM-Newton data for HD 81809 
and carefully determine the presence and signifi- 
cance of any periodicities in its lightcurve. Figure 
[2]displays fitted APEC flux values as well as a LSP. 
Indeed, the periodic behavior has continued with 
what appears to be an ^8 year period. 

The LSP produces a local minimum at the ex- 
pected period surrounded by two peaks at a pe- 



5 



le-14 



Star #7, P = 5.06, x 2 = 1-3 




2000 



2002 



2004 2006 2008 

Observation Date 



2010 



2012 




Fig. 1. — Top: Lightcurve for #7. The dashed line is the best-fit sinusoid (of period P m , the most likely 
peak-period according to MC analysis) through the data. Bottom Left: A cumulative histogram of FAP 
values from the N = 1000 Monte Carlo-sampled lightcurves from #7. If the cumulative histogram of FAP 
values falls within the square in the upper left corner, the lightcurve is periodic by the standards employed in 
this paper. Bottom Right: A non-cumulative histogram of peak periods (inverse of the Lomb-Scargle peak 
frequency) detected in the MC-sampled lightcurves of #7. The vertical dashed line indicates the maximum 
period investigated by the Lomb-Scargle Periodogram, 1.2 x B where B is the baseline for usable observations 
of #7. 
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Fig. 2. — Top: A plot of fitted APEC flux for 
HD 8f809, showing a strong ~ 8 year periodic- 
ity. For better clarity, MOS2 data is not shown 
here but agrees with MOSf data within la. Bot- 
tom: LSP of empirical flux for HD 8f809. An 
S-index cycle p eriod of 8.2 years was found by 
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Baliunas et al.l (|1995l) . The apparent disagree- 
ment between the LSP peak frequency for the X- 
ray lightcurve of HD 81809 and the 8.2 year S- 
index period is likely due to semi-regular observa- 
tion timing. 



riod of 6.2 years and 12 years. The MC sim- 
ulations confirm that the 6 yr peak has a FAP 
of <5% in 100% of the MC steps. The MC sim- 
ulations also show that the shift in period com- 
pared to the S-index observations is an artifact 
caused by infrequent time-sampling of the XMM- 
Newton observations c o mpar ed to the observa- 
tions of iBaliunas et ah ( 19951 ). Similar behavior 
can be found in Fourier analysis of X-ray data 
for 61 Cyg, where the double peak in the Fourier 
power close to the known S-index period is re- 
vealed to be an artifact of infrequent and semi- 
regular observation times: A similar phenomenon 
occurs in the Fourier analysis of S-index data when 
sampled at the same time s as X-ray observations 
([Hempelmann et al. I l2006h . 



3.6. Detection capabilities 

Despite the inaccuracy of the exact cycle pe- 
riod, analysis of HD 81809 shows that our method 
for determining the presence of stellar cycles is sta- 
tistically sound. No significant periodic behavior 
was found in the lightcurves of any of the 9 stars 
in our sample. 

To determine which cycles our method is ca- 
pable of detecting, we use MC simulations again, 
assuming that a source's X-ray flux over time is: 



f\t) — F aV g -\- Famp sin 



/ 2nt 



(1) 



and that the observations are subject to Gaussian 
uncertainty with standard deviation, a. Thus, for 
a given period P and cycle amplitude relative to 
the average flux {F amp / F aV g) — s, a MC simula- 
tion creates T randomly spaced data points over 
a span of S years, produces a LSP and calculates 
the FAP of the best fit sine function. N sam- 
pled lightcurves are created for each (P, s) value. 
We then calculate the fraction of N = 1000 MC 
lightcurves that produce a FAP less than 0.3. Sev- 
eral contour plots are shown in Figure [3J 

Contour plots like those in Figure [3] are cal- 
culated via a similar MC algorithm for each indi- 
vidual source. Instead of sampling N observations 
randomly distributed over the baseline with a con- 
stant uncertainty, the observation dates and flux 
uncertainties for the particular star of interest are 
used to generate sample lightcurves in the MC al- 
gorithm. 
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4. Discussion 
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Fig. 3. — Contour plots of the probability of ob- 
serving an underlying cycle with a Lomb-Scargle 
false alarm probability less than 0.3 for a given 
Fampl F avg and cycle period. This is calculated 
assuming that observations are randomly dis- 
tributed over a baseline of 11 years with Gaussian 
flux uncertainties of a x F avg for each observa- 
tion, top: 20 observations, small flux uncertainties 
(a = 0.15); bottom: 10 observations, a = 0.3. 



The time span covered by the observations dif- 
fers for each analyzed field and ranges between 6 
and 11 years. Most sources have more than 30 
exposures with a usable time longer than 10 ks. 
However, several exposures are often taken within 
a few days and thus do not contribute more infor- 
mation to the periodogram. Typically, about 10 
datapoints with sufficient spacing exist. The sta- 
tistical errors depend on the source flux and the 
exposure time. Source flux varies from source to 
source, and exposure time varies even between ex- 
posures. Here, a typical value of relative error on 
the flux is about 25%. 

We now seek to use our results to put an upper 
limit on the prevalence of coronal X-ray cycles in 
stars of spectral type F through M that have 5-13 
year periods. We calculate P(FAP < 0.3) > 0.68 
contours for each source and use these to create a 
binary function over the P-s domain that is equal 
to 1 when P(FAP < 0.3) > 0.68 and otherwise. 
The binary functions for each star are summed 
to generate a map (Z(P,s)) that represents the 
number of cycles that could have been detected 
using our method for a cycle of particular strength 
and period. 

We assume for simplicity that the intrinsic rate 
of coronal X-ray cycles is approximately constant 
across spectral types F-M. This assumption may 
or may not be justified: th e S-in dex monitoring 
presented in iBaliunas et al. ( 1995 ) did not inves- 
tigate M stars (except for one for w hich no cycle 
was fo und), however, monitoring bv iMauas et al 



( 20121 ) has detected possible S-index cycles in M- 
dwarfs which fall within our 5-13 year domain. 



Baliunas et al.l ([19951 ) do see differences in chromo- 
spheric cycles between F, G and K type stars, such 
that e.g. almost all K-type stars with a low S- 
index have a cycle, whereas F-type stars with low 
S-index typically have flat activity curves. How- 
ever, the weakness of the chromosphere and the 
low contrast of the photosphere make the detec- 
tion of cycles on F-type stars more difficult, so 
that it is not clear if the cycle frequency actually 
varies with spectral type. 

For a particular kind of cycle (P c , s c ), Z(P C , s c ) 
is then the effective sample size. The proba- 
bility of not detecting a (P c ,s c )-type cycle is 
(1 — R) z ( Pc ' Sa \ where R is the intrinsic rate of 
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(P c , s c )-type coronal X-ray cycles for F-M stars. 
A 2a upper-limit on R would imply that R < R 2a 
where R 2a is defined to be such that the prob- 
ability of not detecting a certain type of cycle, 
(1 - R 2 a) Z{Pc ' Sc \ is 0.05. Solving for R 2a in terms 
of Z, R 2a {P,s) = 1 - {Q.Qbf/ 2 . A plot of Z is 
provided in Figure 01 with implied R 2a values pro- 
vided in parentheses. 

While there is a notable area where Rq is high 
(i.e. where the 2a upper limits on the underlying 
rate of coronal X-ray cycles is weak), the statistics 
from this sample strongly limit the prevalence of 
short term cycles (5-8 year periods) in our sample. 
Since both the sun and HD 81809 have s values 
close to 0.7, we can say with 95% confidence that 
fewer than ^50% of the stars in our sample have 
5-8 year X-ray cycles with amplitudes of s = 0.7. 

In addition, if one makes the assumption that 
coronal X-ray cycles should have strengths and pe- 
riods that are more or less uniformly distributed 
over the (Q < s < 1) and (5 years < P < 13 
years) domain, Figure |4] implies a 2a upper limit 
on the ubiquity of coronal X-ray cycles in stars 
represented by our sample. Under this assump- 
tion, our statistics imply that less than ~ 72% of 
those stars have coronal X-ray cycles. 

All statistical constraints given thus far arc 
valid for a stellar population that is well repre- 
sented by our sample, i.e. the population of stars 
detected in a serendipitous X-ray survey. Since 
we are limited by the flux of sources in the field, 
our sample is biased towards X-ray active stars 
(Malmquist-bias). Only two sources have activity 
levels Lx/Lb i < —4.5 and even the least active 
one is still more active than the Sun by an order 
of magnitude or more. Thus, we cannot draw reli- 
able conclusions about F-M stars in general from 
our data. This requires a volume limited sample 
as e.g. in 



( Micela fc Marinoll2003 ) 



Schmitt fc Liefkel (|2004h 



The lack of detection of activity cycles in our 
sample would therefore support the view that in 
active stars coronal X-ray cycles are not common. 
This was observed in the S-index monitoring by 



Baliunas et al 



_ (|1995l ). In X-rays the abse nce of 
cycles in active stars has been suggested bv ISternl 
(jl998l ) on the bases of more fragmentary data. As 



stars age th eir median luminos i ty decreases (see 



reviews of Rosner et al.l Il985t iFavata fc Micela 



2003, and references therein) and they develop 
coronal cycles with larger amplitudes of variation 



5. Conclusions 

We have identified 9 stellar X-ray sources 
that have frequently been observed in repeat 
XMM-Newton observations over the lifetime of 
the observatory. From single temperature APEC 
model fits to extracted spectra, we have gener- 
ated lightcurves of all stellar sources and tested 
the lightcurves of these stars for the presence of 
underlying periodicities using a Lomb-Scargle Pe- 
riodogram assisted by a Monte Carlo algorithm. 
No periodicities were found. 

Using Monte Carlo simulations, we calculated 
which kinds of cycles should be detectable with our 
methodology for each of the 9 stars in our sample. 
We term the underlying rate of coronal X-ray cy- 
cles in F-M stars that have 5-13 year periods as 
R. We assume that R is constant across spectral 
types F through M and that the short-term X-ray 
variability of our stellar sources is negligible com- 
pared to any long-term periodic behavior. Using 
these assumptions, we calculate 2a upper limits 
on R for each type of cycle investigated in this 
study (cycles of any strength that have 5-13 year 
periods) . Our sample is biased towards stars with 
higher levels of X-ray activity; the least active star 
is at least an order of magnitude more X-ray active 
than our sun. If one makes the assumption that 
X-ray active solar-type stars have cycles with pe- 
riods and strengths that are uniformly distributed 
over the (0 < s < 1) and (5 years < P < 13 
years) domain, we can conclude with 95% confi- 
dence that less than 72% of the active main se- 
quence stars have 5-13 year coronal X-ray cycles 
of any strength. We stress that these numbers 
do not apply to the stellar population as a whole, 
which contains a large fraction of stars that have 
significantly lower levels of X-ray activity. 
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and assistance with this project throughout the 
summer of 2011. This work is supported in part 
by the National Science Foundation Research Ex- 
periences for Undergraduates (REU) and Depart- 
ment of Defense Awards to Stimulate and Support 
Undergraduate Research Experiences (ASSURE) 
programs under Grant no. 0754568 and by the 
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Fig. 4. — Contour plot of cycles that could be de- 
tected using our methodology. The z-axis repre- 
sents the number of stars in our sample for which 
a certain type of cycle should have been detected, 
based on MC simulations. 2a upper limits on the 
intrinsic rate (R) of coronal X-ray cycles in F-M 
stars, assuming that R is relatively constant across 
these spectral types, are provided in parentheses. 
The relative amplitude and period of X-ray cycles 
for both the Sun and HD 81809 are denoted by 
solar and star symbols respectively The plot has 
been slightly smoothed for display. 



Smithsonian Institution. We thank an anonymous 
referee for his help in improving the paper. 
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